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1 Introduction 



The aim of the "fastNLO" project is to make the inclusion of jet data into global 
fits of parton density functions (PDFs) feasible. Due to the prohibitive computing 
time required for the jet cross sections using standard calculation techniques, jet data 
have either been omitted in these fits completely or they were included using a simple 
approximation. The fastNLO project implements a method that offers exact and 
very fast pQCD calculations for a large number of jet data sets allowing to take full 
advantage of their direct sensitivity to the gluon density in the proton in future PDF 
fits. This includes Tevatron jet data beyond the inclusive jet cross section and also 
HERA jet data which have been used to determine the proton's gluon density [H 13 
El E] , but which are ignored in current PDF fits El El • 



2 Concept 

2.1 Cross Sections in Perturbative QCD 

Perturbative QCD predictions for observables in hadron-induced processes depend on 
the strong coupling constant and on the PDFs of the hadron(s). Any cross section 
in hadron-hadron collisions can be written as the convolution of the strong coupling 
constant in order n, the perturbative coefficient c„^j for the partonic subprocess 
and the corresponding linear combination of PDFs from the two hadrons Fi which is 
a function of the fractional hadron momenta Xa^b carried by the partons 

a{Hr,llf) = ^ Cn,i{Xa, Xfo, /i,-, /i/) ® [a"(/ir) " Fi{Xa, Xb,Hf)] . (l) 
n,i 

The PDFs and also depend on the factorization and the renormalization scales fif^r, 
respectively, as does the perturbative prediction for the cross section in finite order n. 
An iterative PDF fitting procedure using exact NLO calculations for jet data, based 
on Monte-Carlo integrations of (^, is too time-consuming. Only an approximation 
of (HI) is, therefore, currently being used in global PDF fits. 



2.2 A Simple Approach 

The "/c-factor approximation" as used in [HI |7| parameterizes higher-order corrections 
for each bin of the observable by a factor k = ^^'^ = — computed from 

O'LO 0"(2) 

the contributions with n = 2 (cr(2)) and n = 3 {<J(3)) for a fixed PDF, averaged over 
all subprocesses i. In the iterative fitting procedure only the LO cross section is 
computed and multiplied with k to obtain an estimate of the NLO cross section. This 
procedure does not take into account that different partonic subprocesses can have 
largely different higher-order corrections. Fig. ^ shows that the fc-factors for quark- 
only and gluon-only induced subprocesses can differ by more than ±20% from the 
average. The is therefore minimized under an incorrect assumption of the true 
PDF dependence of the cross section. Further limitations of this approach are: 
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Figure 1: The fc-factor for the inclusive pp jet cross section at ^/s = 1.96 TeV as a 
function of pr at different rapidities y for the total cross section (solid line) and for 
different partonic subprocesses: gluon-gluon (dashed), gluon-quark (dotted) and the 
sum of all quark and/or anti-quark induced subprocesses (dashed-dotted). 



• Even the LO Monte-Carlo integration of (Q) is a trade-off between speed and 
precision. With finite statistical errors, however, theory predictions are not 
ideally smooth functions of the fit parameters. This contributes to numerical 
noise in the calculations |HI distorting the contour during the PDF error 
analysis, especially for fit parameters with small errors. 

• The procedure can only be used for observables for which LO calculations are 
fast. Currently, this prevents the global PDF analyses from using Tevatron dijet 
data and DIS jet data. 

In a time when phenomenology is aiming towards NNLO precision E] , the A;-factor 
approximation is clearly not satisfying concerning both its limitation in precision and 
its restrictions concerning data sets. 

2.3 The fastNLO Solution 

A better solution is implemented in the fastNLO project. The basic idea is to trans- 
form the convolution in (0) into the factorized expression (j^. Many proposals for 
this have been made in the past, originally related to solving the DGLAP parton 
evolution equations P and later to computing of jet cross sections [TIH [TT| IT^ IT^ IT^ . 
The fastNLO method is an extension of the concepts developed for DIS jet produc- 
tion ^01 ^1 which have been applied at HERA to determine the gluon density in the 
proton from DIS jet data jl]. Starting from for the following discussion the renor- 
malization scale is set equal to the factorization scale (/ir,/ = /i)- The extension to 
/ir 7^ /i/ is, however, trivial. The x dependence of the PDFs and the scale dependence 
of a" and the PDFs can be approximated using an interpolation between sets of fixed 
values x^'') and /i^™^ (/c = 1, ■ ■ ■ , fcmax ; m = 1, ■ ■ ■ , mmax) 

<(/i) ■ Fi(Xa,X6,/i) ~ ["=" is true for fcmax,^max,"^max ^ OO] 

■ i^,(4'=\4^/i(")) ■ e^'\x,) ■ e(')(x,) ■ b^^\fi) (2) 

k,l,m 
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Figure 2: Contributions of different partonic subprocesses to the inclusive jet cross 
section at RHIC (left), the Tevatron (middle) and the LHC (right) as a function 
of pt and xt = 2pT/y/s. The subprocess gq jets has been separated into the 
contributions (2) and (3) where either the quark- or the gluon momentum fraction is 
larger. 

where e'^^'^\x) and b^"^\n) are interpolation functions for the x and the fi dependence, 
respectively. All information of the perturbatively calculable piece (including phase 
space restrictions, jet definition, etc. but excluding and the PDFs) is fully contained 
in the quantity 

5-n,i,fc,/,m(/i) = Cn,i{Xa, Xb, IJ.) ® [e^''\xa) " e^^\xb) ■ ft^^H/^)] • (3) 

In the final prediction for the cross section the convolution in is then reduced to 
a simple product 

a(^) ^ J2 a„,,fc,,,4/i) ■ <(/x("^))-F,(xi^-Ui^/x('")). (4) 

n,i,k,l,m 

The time-consuming step involving the calculation of the universal (PDF and 
independent) a is therefore factorized and needs to be done only once. Any further 
calculation of the pQCD prediction for arbitrary PDFs and as values can later be 
done very fast by computing the simple sum of products in While the extension 
of the method from one initial-state hadron ^3] to two hadrons was conceptually 
trivial, the case of two hadrons requires additional efforts to improve the efficiency 
and precision of the interpolation. Both, the efficiency and the precision, are directly 
related to the choices of the points x*-^''-*, /i*-™-' and the interpolation functions e(x), 
b{fi). The implementation in fastNLO achieves a precision of better than 0.1% for 
fcmax, ^max = 10 and mmax < 4. Computation times for cross sections in fastNLO are 
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Figure 3: Comparison of PDF uncertainties for the inclusive jet cross section at RHIC 
(left), the Tevatron (middle) and the LHC (right). The uncertainty band is obtained 
for the CTEQ6.1M parton density functions and the results are shown as a function 
of pt and xt = / ^/s. 

roughly 40-200 /xs per order ag (depending on mmax)- Further details are given in 
Ref dS]. 

The a in (jni) are computed using NLOJET++ jTHllIZI. A unique feature in fastNLO 
is the inclusion of the 0{a'f) threshold correction terms to the inclusive jet cross 
section [TH], a first step towards a full NNLO calculation. 

3 Results 

Calculations by fastNLO are available at http://hepforge.cedar.ac.uk/fastnlo 

for a large set of (published, ongoing, or planned) jet cross section measurements 
at HERA, RHIC, the Tevatron, and the LHC (either online or as computer code 
for inclusion in PDF fits). Some fastNLO results for the inclusive jet cross section 
in different reactions are shown in this section. The contributions from different 
partonic subprocesses to the central inclusive jet cross section are compared in Fig. |21 
for different colliders: For pp collisions at RHIC and the LHC, and for pp scattering 
at Tevatron Run II energies. It is seen that the quark-induced subprocesses are 
dominated by the valence quarks: In proton-proton collisions (RHIC, LHC) the quark- 
quark subprocesses (4,5) give much larger contributions than the quark-antiquark 
subprocesses (6,7) while exactly the opposite is true for proton-antiproton collisions 
at the Tevatron. The contribution from gluon-induced subprocesses is significant at all 
colliders over the whole pt ranges. It is interesting to note that at fixed xt = IptI \fs 
the gluon contributions are largest at RHIC. Here, the jet cross section at xt = 0.5 
still receives 55% contributions from gluon-induced subprocesses, as compared to only 
35% at the Tevatron or 38% at the LHC. As shown in Fig.|21 this results in much larger 
PDF uncertainties for the high xt inclusive jet cross section at RHIC, as compared to 
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Figure 4: An overview of data over theory ratios for inclusive jet cross sections, 
measured in different processes at different center-of-mass energies. Tfie data are 
compared to calculations obtained by fastNLO in NLO precision (for DIS data) and 
including 0{aj) threshold corrections (for pp data). The inner error bars represent 
the statistical errors and the outer error bars correspond to the quadratic sum of 
all experimental uncertainties. In all cases the perturbative predictions have been 
corrected for non-perturbative effects. 



the Tevatron and the LHC for which PDF uncertainties are roughly of the same size 
(at the same xt)- This indicates that the PDF sensitivity at the same xt is about 
the same at the Tevatron and at the LHC, while it is much higher at RHIC. 

An overview over published measurements of the inclusive jet cross section in 
different reactions and at different center-of-mass energies is given in Fig. EJ The 
results are shown as ratios of data over theory. The theory calculations include the best 
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Figure 5: The phase space in x and px covered by the data sets shown in the previous 
figure. 

available perturbative predictions (NLO for DIS data and NLO + 0{a^) threshold 
corrections for pp data) which have been corrected for non-perturbative effects. Over 
the whole phase space of 8 < < 700 GeV jet data in DIS and pp collisions are well- 
described by the theory predictions using CTEQ6.1M PDFs [7j. The phase space in 
X and pt covered by these measurements is shown in Fig. demonstrating what can 
be gained by using fastNLO to include these data sets in future PDF fits. A first 
study using fastNLO on the future potential of LHC jet data has been published in 
Ref. [in]. 
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